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We study the SUSY SO (10) GUT with the symmetric two-zero textures of quark/lepton 
mass matrix which reahze the Georgei-Jarlskog relations in down-type quark and charged 
lepton masses. We show that the important constraints to such framework come from the 
bottom-tau unification and the observed value of sin 2/3, one of the angles in the CKM 
^ . unitarity triangle. We investigate the symmetric two-zero textures assumed at the GUT 

' scale by solving the MSSM renormalization group equations with right-handed neutrino 

(N ■ threshold effects. As a result, the value of tan/3 is constrained to be large enough and the 

representation of Higgs field which couples to the third generation of left- and right-handed 
neutrinos is restricted by the bottom-tau unification condition. The textures with vanishing 
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' 1-2 (and 2-1) element for up-type quarks and vanishing 1-3 (and 3-1) element for down-type 
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quarks are favored by the observation of sin 2/3. 



'"q_i' §1. Introduction 

Q \ One of the most successful predictions in supersymmetric grand unified theories 

' (GUTs) is a gauge coupling unification at the GUT scale. On the other hand, 

any sufficient achievement of matter (Yukawa coupling) unification has not been 
completed yet and it is one of the most important tasks to explain the observed 
quark/lepton masses and mixings based on the GUT framework . 



5^ , Here we would like to notice that the well-known Georgi-Jarlskog (GJ) relation^ 

me = — , m^ = 3ms, rur = rub , 

can provide us with a guiding principle to construct a framework of matter uni- 
fication. The above GJ relations at the GUT scale have been known to explain 
the observed values of down quark and charged lepton masses at low energy quite 
successfully. Especially the first and second relations are quite interesting to re- 
produce the observed down and strange quark masses as well as yu and e lepton 
masses in a unified way. Furthermore we notice that the ratio \/ md/ {rug + fUd)., 
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which is 0.224 it 0.004 experimentally, is almost equal to the observed mixing value, 
iKtslexp = 0.221 — 0.227. This may be a strong indication of zero structure for 1-2 
entry of the down-type quark mass matrix M^, since the contribution of the mixing 
angle from the up- type quark mass matrix M„, which shows more hierarchical struc- 
ture, may be generally expected to be very small. We shall see this coincidence later. 
Thus a natural setup for realizing the GJ relation is to introduce the zero texture of 
symmetric mass matrix with some vanishing entries. Such textures of quark mass 
matrix have been explored by many authors, for example, see 

refs.130 The coeffi- 
cient —3 in the GJ relations is originated from Clebsch-Gordan (CG) coefficients. 

On the other hand, the third relation in eq. (jl-lj) . called bottom-tau unification, 
should be addressed in more detail. Although the predicted value of mb/nir at the 
low energy scale is qualitatively consistent with the experimental value, it depends 
quite strongly on the 3rd generation of right handed neutrino mass Mji^, as well as 
the strong coupling constant 03.^' Actually, it predicts a bit larger value than exper- 
imental values if Mji^ is well below the GUT scale. Furthermore, the recent results of 
the B factory, combined with the great advance in lattice QCD calculations, provide 
us with much more precisely determined observable parameters in the CKM matrix. 
In particular the experimental value of sin 2/?, which is an angle in the CKM unitary 
triangle, has been precisely determined in the last few years. This constrains the 
type of mass matrix textures very strictly. It is shown that the important constraints 
on fermion mass matrix textures come from the measurements of sin 2f^ and that 
the texture with zero 1-3 entry cannot reproduce the experimental data. 

In the previous papers,!^"^ we proposed a symmetric two-zero quark/lepton 
mass texture, which realizes the GJ relations and can successfully reproduce the 
observed bi-large mixings and mass-squared differences of neutrinos in the SUSY 
SO(IO) GUT. In our framework, the predicted branching ratios of lepton flavor 
violation such as — > 67 are safely lower than the experimental upper bound and 
the observed baryon asymmetry of the universe can be explained through thermal 
leptogenesis scenario.'^ It is known (see for example, ref.l^) that ratio of mb/mr 
at the electroweak scale can be reproduced pretty well if M/J3 is of the GUT scale. 
This implies that the two zero texture of neutrino Dirac mass matrix M^^ in our 
model is preferable for predicting mb/rrir because Mr^ is required to be almost of 
the GUT scale in our model, while the right-handed neutrino masses of the 1st and 
2nd generations, Mri and M/j2, are very small. So, we expect that our model can 
reproduce not only neutrino masses and mixing angles but also the down quark and 
charged lepton masses as well as quark mixing angles. However, our model adopted 
the texture with zero 1-3 entry for both M„ and and it is not clear whether it 
would bring enough non-zero contribution to the 1-3 element to be consistent with 
the observed value of sin 2/3. Therefore, the predicted values of sin 2/3 and mb/rrir 
give a crucial information as to whether such matter unification scenarios with the 
GJ relations can be realized in nature, after taking account of the renormalization 
group equations (RGEs) running effect within the framework of unified bottom and 
tau Yukawa couplings. Furthermore, we think it is necessary to find a nice framework 
of zero texture which is consistent with the above current experimental data and at 
the same time which keeps the good GJ relations. 
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The aim of this paper is to investigate the phenomenologically allowed symmetric 
two-zero quark/lepton mass matrix textm'es in the SUSY SO(IO) GUT framework 
with the current experimental bounds of sin 2/3 and mi^/mr. We make numerical 
calculation of the observed values, especially mh/rrir and sin 2/3, by solving the RGEs 
with right-handed neutrino threshold effects. We also show the constraints to Mpts 
and tan /3. *•* 

This paper is organized as follows. In the next section, we present the textures 
of quark/lepton mass matrix which we adopt in this paper. In the third section, 
we show the results of numerical calculations. The final section is devoted to the 
discussion of this paper. 



§2. Symmetric two-zero textures 



We consider the SUSY SO(IO) GUT with the 10 and 126 representation of 
Higgs multiplets. The SO(IO) gauge group is assumed to be broken down to the 
standard model gauge group through the Pati-Salam symmetry. In the above setup, 
the up and down quark mass matrices : M„, M^, charged lepton and neutrino Dirac 
mass matrices : Mg, M^,^ are symmetric ones. Many authors have investigated 
the two-zero textures in the above setup.'^'^ We have assumed that either the 
10 or 126 Higgs multiplet dominantly couples to the fermions in each generations. 
Consequently, we obtain the relation among mass matrices : Mu = M,^^ and = 
Me except for the CG coefficients 1 or —3, which can appear in some elements in 
Me and M^^. The up quark and neutrino Dirac mass matrices are taken to be the 
following textures **^ : 



Ul M„ = I Gu bu Cu I and M^^ = ( x^au X22^« X23Cu | ,(2-1 





XuCtu 






X22bu 


X23Cu 




X23Cu 


X33du 






XnO-u 




X22bu 


X23Cu 


xnttu 


X23Cu 


X33du 



U2 M„ = I bu Cu I and M^^ = \ X22bu X23Cu \ , (2-2) 

where au,bu,Cu and du are complex numbers, and CG coefficients Xij can be taken 
as 1 or —3. We denote the above assumptions of textures as Ul and U2 hereafter. 
We have adopted the texture Ul in refs.'^'^^''^ In this paper, we take xi2 = xis = 
X22 = 2^23 = 1, because the values of Xij except for 2:33 do not affect to the RGEs 
running of the observables in quark and charged lepton sector. The constraint on 
X33 from the bottom-tau unification will be discussed in the next section by showing 
numerical results. 



*■* We write here tan /3 for the notation of the parameter expressing the ratio of VEVs of up and 
down Higgs fields to discriminate it from sin 2/3. 

**-* In addition to the texture Ul which we have adopted in our previous papers,!^^^ we here 
also introduce texture U2. The latter texture is obtained by exchanging the 2nd and 3rd generation 
suffices, and it is expected not to make so much difference at least to the calculated neutrino masses 
and mixings because in the neutrino sector the 2nd and 3rd generations are almost maximally. 
Nevertheless, we guess that it does bring much difference to quark sector. 
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The down quark and charged lepton mass matrices are given by 

ad \ (ad 
Dl Md=\ ad bd Cd and Me = -36^ q | , (2-3) 

Cd dd / \ Cd dd 

where a^, bd, Cd and dd are also complex numbers. The CG coefficient in 2-2 element 
of Me is the crucial ingredient to realize the GJ relations. In addition to the above 
texture, we also consider the following three-zero texture for the down quark and 
the charged lepton mass matrices : 

ad \ I ad 

Dl' Md= \ ad bd 1 and Me = | -36^ | , (2-4) 

dd 

which is a special case of Dl, imposing further zero to the 2-3 element *-*. The texture 
Dl' were originally proposed by Georgei and Jarlsko^ and also investigated in ref.l^ 
The right-handed neutrino mass matrix is taken to be 

aR 

Mr= \ aR I , (2-5) 

dR 

where we assume \aR\ <C \dR\ = Mr^ ~ Mqut where Mr^ is the third generation 
of right-handed neutrino mass. There are two characteristic features of this form. 
First, it reproduces the observed bi-large mixings of neutrinos .'^'^ Second, it pre- 
dicts almost equal Majorana masses for the 1st and 2nd generations of right-handed 
neutrinos. The former is important to derive m-r/rrib ratio, while in the latter, 
degenerate masses guarantee the baryon asymmetry of the universe through lepto- 
genesi^' in taking account of the resonance effects.^ The left-handed Majorana 
neutrino mass matrix is obtained by the seesaw mechanism : Mi, = — MJ^M^^'^Mj^^ . 



§3. Numerical analysis 



In this section, we show the predicted values of the observable parameters in the 
CKM matrix : VckM) the quark and charged lepton masses at the electroweak scale 
by solving RGEs with the right-handed neutrino threshold effects. 

The procedure that we adopt here for generating scattered plots of Fig. |21 Fig. |31 
Fig. m and Fig. El is explained as follows. First, we generate random numbers for 
au,d,bu^d,Cu,d:du^d ill the mass matrices. Then, after solving the RGEs for Yukawa 
couplings with the generated values of a^^dibu^, c^^d-, du^d at the GUT scale, we plot 
the observable parameters if the up quark and charged lepton masses fall into the 

*■* It might be possible to introduce the case, in which we impose further zeros to the matrix 
elements of M„. However, we have already recognized that it does not reproduce the neutrino bi- 
large mixing angles any more and the minimum parameter set should not be less than 4 for M,y„ . 
Thus the textures, Ul and U2 with Dl and Dl', are natural extension of our previous modeP*""! 
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Fig. 1. (a): The as and tan/3 dependence of mt and are shown. The heaviest right-handed 
neutrino mass is taken to have the range of Mrz = 10^ ~ 10^ GeV. (b) : The as and Mr3 
dependence of rub and m-r are shown. The tan/3 is taken to have the range of tan/3 = 3 ~ 50. 
The minimal and maximal value of as are taken according to ea. (l3-2ll in both (a) and (b), 
respectively. The lighter right-handed neutrino masses are fixed as Mm — Adm = 10^" GeV in 
these figures. 



following experimentally allowed range : 

rUuiMz) = 0.8 - 3.0 MeV,mc{Mz) = 500 - 800 MeV,mj(Mz) = 170 - 180 GeV, 
me{Mz) = 0.487 MeV, m^{Mz) = 103 MeV, mriMz) = 1.75 GeV, 
iKsl = 0.221 -0.227 , (3-1) 

where the m^, rUc and mt are taken to have maximally wide range,^ the charged 
lepton masses are allowed to have 1% error around the displayed central value.C^ 
We then obtain the predicted values of ma, nig, m^, \Vub\, \Vcb\ and sin 2/3. The 
another important observable is the strong coupling constant which is defined as 
"3 = 5-3 /4vr : 

asiMz) = 0.1187 ± 0.0020 , (3-2) 

where we use the value in.^^- As we will see later soon, the low-energy value of the 
bottom to tau lepton mass ratio depends on the value of 03. 

3.1. bottom-tau unification 

Here, we discuss the constraint which comes from bottom-tau unification. Before 
showing the numerical results, we present the parameter dependence of the bottom 
to tau mass ratio at the low-energy scale. In the MSSM with heavy right-handed 
neutrinos, the 1-loop RGE for the ratio of bottom and tau masses is approximately 
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Fig. 2. (a): The scattered plots of mt and rrir in the case of Ul-Dl textures for tan/3 — 3 and 50. 
The X33 is fixed as 3:33 = 1 in this case, (b): The case of 1:33 — 1, —3 are shown when tan/3 is fixed 
as tan/3 — 50. In both figures, the right-handed neutrino masses are fixed as Mm = Mr2 = 10^" 
GeV and Mrs = 10^^ GeV. The two horizontal lines in the figures correspond to the allowed 
range of mr. The vertical lines correspond to the allowed range of mt- We find that the case of 
large tan/3, say tan/3 ~ 50 and a;33 = I is favored by bottom-tau unification comparing with the 
experimentally allowed range for bottom quark mass : 2.8 < mj, < 3.1 GeV from these figures. 



given by 

d ( nib \ If rub 



dt J IGvr^ \m 



iyu)h - iyu)h + 3 ((yd)33 - (ye)i3) - ( - ^al 



(3-3) 



where 31^3 are the gauge couphngs of the standard model gauge group, yu,d,i^,e are 
the Yukawa couphngs for the up and down quarks, neutrinos and charged leptons, 
respectively, and t = \nfi. The top quark Yukawa coupling {yu)33 tends to cancel 
the contribution of the 53 term on the right hand side in ea. ()3-3() . This cancellation 
is important to reproduce the phenomenologically allowed low-energy value of the 
bottom to tau mass ratio. Since the tau neutrino Yukawa coupling (2/1^)33 cancels the 
top quark one, the bottom quark mass becomes smaller if the absolute value of the 
CG coefficient X33 in 3-3 element decreases. Moreover, the lower values of tan f3 and 
Mji3 which determines a decoupling scale of the third generation of right-handed 
neutrino are also disfavored for the same reason. The explicit example of plots in 
m-b — rrir plane are shown in Fig. ^for the different values of 03, Mpts and tan/3, 
respectively. In these figures, the RGEs are numerically solved. The 03 and tan /? 
dependence of rub and nir are shown in Fig.^a). The heaviest right-handed neutrino 
mass is taken to have the range of Mj^ = 10^^ ~ 10^^ GeV in Fig. ^a). The as 
and Mfi3 dependence of rub and nir are shown in Fig. E^b) . The tan /3 is taken to 
have the range of tan f3 = 3 ^ 50. The minimal and maximal values of are taken 
according to eq. H3-2() in Fig. ^a) and^b), respectively. It is shown that both mb 
and nir are proportional to tan/3. Furthermore, nib ("ir) is also proportional to 
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as (M/js), respectively. We can confirm the above arguments by the figures of the 
results of numerical calculations. 

Fig. shows the scattered plots in — m,- plane in the case of Ul-Dl for 
the different values of tan /3 and X33 in taking account of the constraints shown in 
eqs. H3-l|) and ()3-2|) except for rrir- In this analysis, the right-handed neutrino masses 
are fixed as Mm = Mr2 = 10^° GeV and Mm = 10^^ GeV, respectively. In the 
Fig. Efa), we take 3:33 = 1. We find that large tan/3 is favored by the bottom-tau 
unification with allowed range of bottom quark mass : 2.8 < < 3.1 GeV. The 
Fig. Hfb) shows the X33 dependence of bottom quark mass. We easily find that the 
allowed values of bottom quark mass are obtained in the case of X33 = 1. Therefore, 
we will consider only the case of 

X33 = 1 and tan /3 = 50 , (3-4) 

hereafter in this paper. 

Our next task is to show the difference of the down and strange quark mass 
predictions in the textures of mass matrices. Since the low energy value of the bottom 
to tau mass ratio is not sensitive to the textures of the up quark and neutrino Dirac 
mass matrices, it is enough to show the difference in the case of Dl and Dl'. For the 
case of Ul-Dl textures with X33 = 1 and tan /5 = 50, the down and strange quark 
masses are predicted to be 

md(Mz) = 3.7 - 5.0 MeV , ms{Mz) = 55 - 74 MeV , (3-5) 

where the constraints of eas. H3-l|) and (|3-2() are considered, respectively. Almost the 
same values are predicted in the case of U2-D1. These values are nicely located within 
the experimental allowed region : md{Mz) = 2.6 — 5.2 MeV and ms{Mz) = 52 — 85 
MeV.I^'^ For the case of Ul-Dl' textures with X33 = 1 and tan/3 = 50, the down 
and strange quark masses are predicted to be 

md{Mz) = 3.2 - 3.4 MeV , ms{Mz) = 81-85 MeV . (3-6) 

Almost the same values are predicted in the case of U2-D1'. These values are more 
strongly restricted rather than the case of Dl. We recognize that these two textures 
of down quark mass matrix, M^^ can reproduce the observed down and strange quark 
masses. 

3.2. The CKM matrix and unitarity triangle 

The recent measurements of the CP violating decay of B mesons into charmo- 
niums provide us with the precise value of sin 2/3.^^' 

sin 2/3 = 0.685 ±0.032 , (3-7) 

where (3 is one of the angles in the CKM unitarity triangle. The angle (3 is the 
most precisely measured quantity among the angles in the CKM unitarity triangle 
at present. The angle (5 is given in terms of the elements in CKM matrix by 
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Fig. 3. (a): Predicted values of |Vc6| are plotted as a function of \Vub\ for the case of Ul-Dl. (b): 
Predicted values of sin 2/3 are plotted as a function of |Ki6| for the case of Ul-Dl. In both figures, 
the constraints shown in eas. ll3-ll l and 13-211 are considered. The relevant parameters are fixed 
as tan^ = 50, 3:33 = 1, Mm = Mr2 = 10^° GeV and Mr3 = lO^'^ GeV, respectively. The two 
horizontal lines in the figures correspond to the allowed range of \Vcb\ and sin 2/3, respectively. 
The vertical lines correspond to the allowed range of |Vui,|. 



We also take the experimental allowed range of the CKM matrix elements^' : 

\Vub\ = 0.0029 - 0.0045 , \Vcb\ = 0.039 - 0.044 . (3-9) 

It is obvious that the predictions of sin 2(3 and the CKM matrix elements are closely 
related each other. 

Let us show the difference between the cases of Ul and U2 at first. Before 
showing the numerical results, we estimate the absolute values of the CKM matrix 
elements and sin 2/3 by following the discussion presented in ref.l^ For the case of 
Ul-Dl, we obtain the following approximated relations : 

Ks -W — , Tf^ - \/ ' /3-arg 1-W e , (3-10) 

where the CP violating phase (p is the combination of ones of the quark mass matrices. 
These approximate forms are given at the GUT scale. One expects almost constant 
ratio of \Vub\ to \Vcb\ from ea. ()3-lU|) . For the case of U2-D1, we obtain the following 
relations : 




rriu ^ f, Imurus 1 j^/ 



\Vus\ ^ \ — , \Vub\ — , /3 ^ arg 1 - ^ , (3-11) 

mt V V '^trrid \Vcb\ 



where the CP violating phase (j)' is also a combination of ones of the quark mass 
matrices. Comparing eq. H3-1U() with eq. H3-ll() . one finds that the angle /3 in the 



Can symmetric texture reproduce unitarity triangle and rrHf/m-r- ? 



9 



0.003 



(a) 

0.0034 0.0038 



0.0042 



(b) 

0.041 0.042 



0.043 



1. - 
0.9 
0.8 
0.7 - 
0.6 

sin2/3 0.5 
0.4 - 
0.3 
0.2 
0.1 



^^^^ 



• • • •» 



.••••V 



4. 



0.9 
0.8 
0.7 
0.6 



•J,r J- 0.5 sin2yS 0.5 



0.4 
0.3 
0.2 
0.1 





0.003 



0.0034 0.0038 
|V„bl 



0.0042 



0.04 



0.041 0.042 
IVcbl 



0.043 



0.044 



Fig. 4. (a); Predicted values of sin 2/3 are plotted as a function of \Vub\ for the case of U2-D1. (b): 
Predicted values of sin 2/3 are plotted as a function of \Vcb\ for the case of U2-D1. In both figures, 
the constraints shown in eas. ll3-lll . 13-21 1 and 13-911 are considered. The relevant parameters are 
fixed as tan/3 = 50, 133 = 1, Mri = Mr2 = 10^° GeV and Mm = lO^'^ GeV, respectively. The 
two horizontal lines in the figures correspond to the allowed range of sin 2/3. 



case of U2-D1 is predicted to be larger than the Ul-Dl and may be favored by 
experimental data. 

Fig. Elshows the predicted values of \Vub\-, and sin 2/3 for the case of Ul-Dl 
which are to be consistent with the constraints shown in eqs. H3-l|) and ()3-2() . In all 
analyses of Fig. |21 Fig. |1] and Fig. |SJ the relevant parameters are fixed as tan /3 = 50 
and 2:33 = 1, Mri = Mr2 = 10^° GeV and Mm = 10^^ GeV, respectively The \ Vuh\ 
is predicted to be near the lower bound in the case of Ul-Dl, and the predicted value 
of \Vch\ is proportional to \Vub\- One can also easily find that the predicted value of 
sin 2/3 does not reach the experimental lower bound. Thus, we should exclude the 
case of Ul due to the observed value of sin 2/3. Therefore, we will investigate only 
the case of U2 in more quantitative way in the following. 

Fig. |1] shows the predicted values of the same parameter sets for the case of 
U2-D1 in taking account of the constraints in eq. ()3-9|) in addition to eqs. ()3-l() and 
H3-2|) . The \ Vub\ and \ Vcb\ are found to cover the whole experimentally allowed region 
in this case. The predicted values of sin 2(5 can reach the experimental allowed region 
only for \Vub\ > 0.0036. We can conclude that there exists the allowed region where 
all the experimental constraints including sin 2/3 are satisfied in the case of U2-D1. 

For the case of U2-D1', we find the predicted values of the same parameter sets 
as in the case of U2-D1, which are shown in Fig. 1^1 This case has the much smaller 
allowed region rather than the case of U2-D1. Since the \Vub\ is constrained to have 
very small value, it seems to be difficult to reach the allowed region of sin 2/3. The 
predictions for the all textures which we take are summarized in the Table 1. We 
can see that the case of U2-D1 is the most favored texture, especially if we compare 
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Fig. 5. (a); Predicted values of sin 2/3 are plotted as a function of \Vub\ for the case of U2-D2. (b): 
Predicted values of sin 2(3 are plotted as a function of \Vcb\ for the case of U2-D2. In both figures, 
the constraints shown in eas. ll3-lll . II3-2I I and 13-911 are considered. The relevant parameters are 
fixed as tan/3 = 50, 133 = 1, Mm = Mr2 = 10^° GeV and Mrs = lO^'^ GeV, respectively. The 
two horizontal lines in the figures correspond to the allowed range of sin 2/3. 





md{Mz) 


ms{Mz) 


mb{Mz) 


\Vub\ 


\Vcb\ 


sin 2/3 


Ul-Dl 


3.7-5.0 


55 - 74 


3.0-3.2 


~ 0.003 


all 


< 0.6 


U2-D1 


3.5-4.8 


55 - 77 


3.0-3.2 


0.0036 < 


all 


< 0.8 


U2-D1' 


3.2 - 3.4 


81-85 


3.0-3.2 


< 0.0034 


all 


~ 0.5 -0.6 



Table I. The table shows the summarized predictions for the different textures in the case of tan /3 = 
50, X33 = 1, Mfli = Mr2 = 10^" GeV and Mrz = lO^'^ GeV, respectively. The "all" means that 
the predicted region covers all the experimentally allowed region. 



the predicted values of \Vub\ and sin 2/? of the other cases. 



§4. Discussion 



We have investigated the symmetric two- zero textures for all quark/lepton mass 
matrices in the SUSY SO(IO) GUT, in which the Georgi-Jarlskog relations are real- 
ized. The low-energy predictions are obtained by solving the renormalization group 
equations with the right-handed neutrino threshold effects. One of the important 
consequence of such framework is that the bottom quark and tau lepton masses are 
unified at the GUT scale. The large tan/3 and the heaviest right-handed neutrino 
mass : near the GUT scale are favored by the bottom-tau unification. The CG 
coefficient in 3-3 element of neutrino Dirac mass matrix : X33 is also constrained to 
be 1 not —3. The predicted value of sin 2/3 strongly depends on the choice of up- 
quark mass matrix textures. As results of the analysis, the case of U2-D1 is favored 
by the measurement of sin 2/3 while the other cases which include Ul-Dl, U1-D2 and 
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U2-D2, are not. However, the allowed region on the \Vub\ — sin 2/3 plane in the case 
of U2-D1 is very limited. 

Finally, we would like to discuss the possible textures and patterns of the CG 
coefficients for the neutrino Dirac mass matrix. These are directly related to neu- 
trino oscillation observables. The threshold effect of the lightest and the next-lightest 
right-handed neutrino masses : Mri and Mr2 in the RGEs running become much 
important to predict neutrino oscillation observables. We can obtain a strong con- 
straint to the parameters in the neutrino Dirac mass matrix and the right-handed 
neutrino one from such analysis. We can expect that such investigation leads us to a 
complete framework of matter unification. 
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